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Abstract 

We demonstrate that the long distance contribution to 5^ — ^ p^7 may 
have a large effect on the decay rate and through it on the CP violating partial 
rate asymmetry. The change in the asymmetry can be between a factor of 
0.4 and 1.4 with respect to the effect of the short distance contribution alone, 
depending on the values of the CKM parameters p and rj . 



1. Introduction 

Interest in rare B decays lies mainly in their potential role as preci- 
sion tests of the Standard Model (SM) . Measurement of the relevant decay 
rates will provide useful information about the Cabibbo-Kobayashi-Maskawa 
(CKM) DH matrix elements. In particular jV^^I is a parameter of crucial im- 
portance in the SM and it is still very poorly known. It is unrealistic to 
expect a measurement of \Vtd\ from t-quark decays in the near future. One 
would then hke to find alternative processes to determine \Vtd\, e.g. pro- 
cesses that are dominated by virtual t-quarks. The present constraints on 
\Vtd\ are derived from the experimental value of the — mixing param- 
eter Xd = AM/r. The mass difference AM is calculated from the usual box 
diagrams. Even though the value of the t quark mass is known 0], this cal- 
culation is still hindered by the uncertainty in the decay constant /b- \Vtd\ 
can be determined precisely from the decay tt+z/z/, which is however 

very rare. 

As is well known, short distance ( SD or penguin) contributions play a 
leading role in b —>■ s{d)'y decays, whose amplitudes are proportional to Vts 
and Vtd respectively. 

Since the decays B K*^ have already been observed it is useful to 
understand what we may learn about CKM elements through a measurement 
of -B — >• p7 decays. 

Recent investigation of the long distance (LD) quark level specta- 
tor contributions to these decays via vector meson dominance ( VMD ) 
{b s{d)V — i> s{d)'y) have shown that they can be neglected. LD contribu- 
tions to the decay b s((i)7 via cc states are only a few percent relative to 
the penguin contribution because of a strong dependence of the ip decay 
constants ( < 0\c'j'^c\iIj >= ie':^g^{q^) ) g^ijn^) » g^{0) . The LD contribu- 
tions to the decay b ^ d'-y via p and a; are only about 3% ■ (y*ij/Vtd) relative 
to the penguin contributions 0. 

Another source of LD contributions to b ^ d'j could come from penguins 
when we take into account the difference between u and c quark masses. How- 
ever Soares has claimed that this contribution is only about 7% ■ (V^V^d) 
of the penguin, so we will neglect it. On the other hand LD contributions to 
the B^ p^'j decays that are not of the spectator type may be significant 

i- 

In this paper we will investigate B^ —> decays, including the W 
exchange diagram (ED) in our analysis (see Fig. la). We will show that the 



1 



ED amplitude is about ±60% ■ V*^^/Vtd of the penguin (SD) amplitude (see 
Fig. lb). The decay rate for — > p^7 can vary by a factor of 0.7 — 2.5 ( with 
respect to the penguin contribution alone ) depending on the CKM matrix 
elements V^f,, Vtd and on the relative sign of the W exchange and penguin 
amplitudes. 

We remark that ED and other non-penguin contributions to B p'j 
were also studied in some detail in refs. ^ 0. However, the approach of 
ref.Q suffers from significant double counting in that it includes ED and 
VMD contributions additively. It is claimed wrongly that ED and VMD 
contributions emanate from two different quark diagrams. The formalism of 
ref.0 basically coincides with ours regarding the ED but differs essentially 
from our treatment of penguin and related contributions. For instance the 
gluons that appear in their lowest order penguin diagrams are absorbed in 
our case in meson wave functions. 

Measurement of CP violation in B meson systems would complement ex- 
isting results in K decays and help elucidate the origin of this phenomenon. 
Observation of the direct CP violation in B decays would eliminate the su- 
perweak model for CP violation. 

The main focus of our paper will be the CP rate asymmetry in the ex- 
clusive decays B^ p^'j: 

_ r(i?+ ^ p+7) - T{B- ^ p-7) 
""^P - r(5+ ^ p+7) + T{B- ^ p-7) ■ ^ ^ 

which was investigated in ref. without taking into account ED contribu- 
tions. 

An important finding of our investigation is that the exchange diagram 
does not contribute significantly to the rate difference of the decays of charge 
conjugate channels B^ (the numerator of eq.(P). As was mentioned above, 
the rate ( i.e. the denominator of eq.(|l|)), can change by a factor 0.7 — 2.5. So 
the CP asymmetry for the decays B^ p^7 will change by the same factor 
as the rates (compared to the SD contribution by itself ). The resulting CP 
asymmetry could be as large as 30% within the allowed range of CKM matrix 
elements. 

Because the ED contribution is negligible for the processes B^ K*^'j 
due to the relevant CKM factors, it can also importantly affect the interesting 
ratio 
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2.W Exchange Diagram contribution to p^7 

It is well known that in some processes ED can be appreciably large i.e. of 
the same order or even larger than the SD (penguin) see Fig. lb, contribution 

0]. The exclusive decays 5^ , can occur via ED as shown in 

Fig. la. In a simple constituent quark model for the B and p mesons ( and 



neglecting terms of order M^/M^)^ the amplitude is given by g, |T0 



Aed[B 7 = ^TT ^eKbKd 

QMBrriu 

Tr[(A-MB) A^^M^^)] (3) 

The QCD correction factor [C2(Mb) + Ci(Mb)/3] ~ 1 has been absorbed, ruu 
= 330 MeV is the constituent mass of the u quark and we used the following 
definition for the decay constants fs and fp{fp^ 216MeV ) 

< 0\uj^d\p- >= le^pMpfp, < 0\u^^^,b\B- >= ifsP^. (4) 

The factor rriu in the denominator of eq. (|^) indicates clearly that the ED 
is a long distance contribution. The effective value for that should be 
used could easily be 30% larger or smaller than the commonly used value 
quoted above. This then constitutes an intrinsic uncertainty in our estimate 
of the ED amplitude ( see however a comment at the end of the paper ). 

Evaluating the trace in (|) we get 



A (U- - \ V2GFfBfpMp 



We denote by pv {p-y) and ey (£7) the momentum and polarization of the 
vector meson (photon), respectively. 



3. Short Distance contribution to B ^ p^y and B K*'y 
Now we consider the SD contribution to the B V'j decay (V = p, K*). 
We will follow refs. |TT[ for the calculation of the exclusive matrix elements 
for B V^7, adding SU (3) flavour breaking effects for the vector meson wave 
functions pl, H. 

The SD Hamiltonian for the decay b —>■ s{d)'j, is 

AC 

H = -—^v,C,0-j, (6) 

where q = s oi d and Vg = VtbV*^, Vd = VibV^^. The magnetic dipole operator 
O7 is (in the limit ^ ^ 0) 
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^7 = jT—^qL^T^^nibbRF,,, (7) 

and we will use 6*7(5 GeV) = —.305 @]. The amplitude for the exclusive 
decay B —>■ Vj is then 

Asd{B ^ V^) = -^v,C7 < V^\0^\B > (8) 

The transition matrix element of the magnetic operator O7 for an on-shell 
photon has the general form (we neglect the vector meson mass) 



< V^\Or\B >= eMBFviOr){te^,ap£%£yp 
In the model of ref.[§] F{0-j) is given by 



FviOj) = — / dy—==(j)v{y)(l)B{ • 10 

4:7!-'^ Jo y/l + y 2 

Here (psip) = exp{ — p^/(2p|n)) is the harmonic oscillator wave function. 
Cb = VStt^^'^Pp^^'^ and Cy = /y/(4-\/3) are the normalization factors of the 
B meson and vector meson wave functions, respectively, (pviv) is the quark 
distribution amplitude in the vector meson ( ps(d) = PvV and Pu = pv{^ — y))- 
For the p and K* mesons we take, respectively ||13 
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<P{y)p = Qy{l -y){l~ 0.85[(2y - 1)^ - -]} (11) 

and 



<t>{y)K' = 6i/(l-y){l + 0.57(2y-l)- 

1.35[(2y - 1)^ - + 0.46[^(2y - If - {2y - 1)]}. (12) 

A.The decay rate for B"^ — >■ p'^'j 

Before presenting our results we would like to note that the relative sign 
of the ED and penguin amplitudes is not known. We will assume that ED 
has no strong phase. Such phase could arise for example from order QCD 
corrections which however are expected to be small because as{mb) ^ 0.2. 
Adding eqs. (^ and (|^) we get 



A{B- ^ p-7) 



V2 



A 



SD 



where 



AsD + Aed fsffMp Vub^ 



AsD ' 6Mim„C7F,(07) V,*, 

The rate for B~ p~7 is then given by 



T{B- ^ p-7) = 2GlMlaC^,\Fp{0-j)\'\Vu\'^ 



where 



AsD + Aed 



A 



SD 



(13) 



(14) 



(15) 



(16) 



As is seen from eq.(^ the ED contribution is proportional to the B meson 
decay constant fs- The B"^ — B^ mixing parameter Xd is proportional to 



q2 

Xd = TB-^\Vtd\'mlMBflBBr]QCDF{ml/M^) 

071 



(17) 
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where rjQCD is a QCD correction factor and 



Fix) 



1 r 3 — 9x \nx 
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fa; -11 



+ 



Here x = m^/M^r and 



m 



pole 



mt 



where i^^t ~ 11 |T^ . 

Substituting fs from eq.(|l^) in eq.(0) we get 



(18) 



(19) 



AsD + A 



ED 



A 



{IT: 



fpM, 



ub 



XdQn'^ 



SD ' mlm^CjFp{Oj) \Vu\V;,\ rBGlmlMBBBr]QCDF{mj/Ml,) 

(20) 

Using the Wolfenstein parametrization for the CKM matrix elements [|l6l 
we obtain 



AsD + Aed _ r. , p- -r]"^ -IT] ^ 



SD 



[(l-p)2+r/2]3/2^ 



where 5^ is given by 



(21) 



9 = T: 



fpMp 



xaQn"^ 



(22) 



For our estimates we use the following input parameters |Tj]: Xd = 0.67 ± 
0.10 , tb = (1.49 ± 0.04)10-^2^, mf'" = 176 GeV § ( = 167 GeV ) , 
riQCD = 0.85 ± 0.05, Bb = I , X = 0.22, A = 0.785 ± 0.093, m„ = 330MeV. 
For the central values of these parameters we get \g\ = 0.6 . 

For Q in the Wolfenstein parametrization of the CKM matrix we have 



AsD + A 



ED 



A 



SD 



{1+9 



p — p^ — if 



9W 



[(l-p)2 + 772]3/2J [(l_p)2+^2]3 



(23) 
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In Fig. 2 we plot curves corresponding to constant values of for g = 0.6 
(Fig. 2a) and g = —0.6 (Fig.2b). Note that if the ED is not taken into account 
{ g = 0) , then f2 = 1 for all p, rj . 

We indicate in the Figs. 2-5 the region in the p, t] plane that is allowed by 
present experimental bounds. The constraints on p and rj are obtained [|14| 
from the following experimental results 



ub 



cb 



0.08 ±0.03 



\€k\ = (2.26 ±0.02) ■ 10" 



(24) 



The bounds on 



circles centered at (0,0' 



give \/~(F^rri^ = -^^(O.OS ± 0.03) which leads to two 



The bounds on ex lead to the hyperbolas [14 



- = 1.5 ■ A^Bk + 6.6 ■ A'^Bk ■ (1 - p) 
V 



(25) 



where we will take A = 0.785 ± 0.093, Bk = 0.8 ± 0.2 |n 
For the branching ratio of the decay B~ — > p~7 we get 



Br(B- 



p-7) = 2TBGlMl,aC',\F,iOr)\'A'\%l-pr + v'Mp,v) 
~ 8.2-10-^[{l- pf + r]^]n{p,r]) ■ 



f_^Y (- 



-Y (^Y 

'/ Vn.785/ 



.-0.305^ V-2.6-10-3^ vo.785^ 
For the ratio R in the Wolfenstein parametrization we get 



(26) 



R 



T{B- ^p-^)+T{B+ ^p+^) 



T{B- ^ 
Fp{Or) 



K* 

2 



-7) ± T{B+ K*+-i) 



Fp{07 



where 



FpjOj) 12 



0.545 ± .025 for pp = 0.5 - 0.65 GeV. 



iv; 



td\ 



Wt 



ts 



n{p,ri) 
(27) 



IFk*(07)I 

As an illustration we plot curves corresponding to constant values of R 
for g = ±0.6 (Fig. 3a) and g = —0.6 (Fig. 3b) in the p,r] plane. 
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5. CP asymmetry in the decay — >■ 

We now consider direct CP violation in the decay B"^ P^l- 
In the spectator approximation the absorptive contribution to the quark 
decay b d'y comes from the real intermediate states b — > uud d'y, 
b — >■ ccd — s> d'y, b gd ^ d'y. The absorptive part of Hamiltonian in the 
spectator approximation is 



AC 1 12 
H,, ^ -^a,-^0r{-C2 ■ K.Kd + 0-12-^2 ■ V^tV;, + -Cs ■ Vt,V;,) (28) 



Here ^2(5 GeV) = 1.096 |[ and 6*8(5 GeV) = -0.185 |T8l. 

As explained earlier, the absorptive part of the ED amplitude is of order 
ttg. Moreover the interference of the ED amplitude with the first term in Hab 
will not contribute to the numerator of eq. (|I|) either. We have also esti- 
mated the contribution stemming from the interference of the ED amplitude 
with the last two terms in eq.(p8D and found that it enhances (reduces) for 
g = —0.6 (^^=+0.6) the numerator of acp by approximately 12% with respect 
to the leading contribution alone. We will therefore ignore it together with 
the order contributions. We conclude that the ED can affect significantly 
only the denominator of Ocp in eq. (|1|). To estimate the numerator of the 
asymmetry in eq.([l|) we follow the formalism of which includes two types 
of contributions. The first results from the interference of the first two terms 
in Hab with the leading penguin amplitude ( see eqs.(|^)-(|^)) while the second 
includes the interference of the absorptive parts in non-spectator diagrams 
(see Fig. 3 in ref. [Q) with the leading penguin amplitude. Since the ab- 
sorptive part comes from long distance effects, a cutoff parameter A for the 
gluon momentum, is introduced to avoid double counting with the contribu- 
tion from the wave functions of the mesons [§. The denominator of eq.(|I]) is 
calculated using eq.(|T|) {T{B- p-7) + T{B+ p+7) ~ 2T{B- p-7)). 

The asymmetry a^p is then given by 



a. 



cp 



- -^A .. 1,.. ^^.^ (29) 



[(l_p)2 + ^2] 2 Crn{p,v) 



where Ka ^ 0.97 for A = 1 GeV and Ka ~ 1.3 for A = 0.2 GeV [|]. 

In Fig. 4 we plot curves of constant acp ( a few representative values) for 
g = +0.6 (Fig. 4a) and g = —0.6 (Fig.4b) for cut-off parameter A = 0.2 GeV 
and a, = 0.2. 



8 



We now estimate the total number of charged B decays required to ob- 
serve this CP asymmetry at the 3a level, N^u, assuming perfect experimental 
detection efficiency 



a%Br{B 7) 
^ 5 ■ 10^ ^i^l^^^ ■ fi(p, r^) ■ (30) 
/1.3n2/0.2x2/-2.6 • 10-3s 2 /0.785n2 

In Fig. 5 we plot curves of constant N^a- for g = +0.6 (Fig.Sa) and g = —0.6 
(Fig.5b). 



Summarizing our main results, Fig. 4 indicates that there can be a sizable 
CP asymmetry in P^l, possibly as large as 30%. Comparison of Figs. 

4a and 4b shows that the long distance ED amplitude plays an important 
role in determining the value of the asymmetry ( a factor 0.4-1.4 with respect 
to the pure penguin contribution), and also underlines the need to theoret- 
ically understand the relative sign of the penguin and ED amplitudes. The 
importance of the annihilation diagram is also apparent from Figs 3a and 3b, 
where the ratio R defined in eq. (0) is shown. The plot of Fig. 5 shows that 
about 10^ — 10^ charged B mesons are needed to be able to observe the CP 
partial rate asymmetry at the 3a level under ideal experimental conditions. 
Such numbers are not out of the question for near future B factories. Finally, 
let us note that the magnitude of the ED contribution as presented here will 
be determined by checking against experiment the predictions of ref. for 
B — > where the ED is supposed to dominate the amplitude, leading to a 
branching ratio of about 4 ■ 10~^ for rriu = 330 MeV. 
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Figure Captions 



Fig. la The W Exchange Diagram (ED) for the decay B ^ p'y 

Fig. lb The electromagnetic penguin diagram for the decay B V'j 

iy = K\p). 

Fig. 2a Curves corresponding to constant VL = 0.75,0.85,1.0,1.5,2.0,2.5 
for g = 0.6. The full lines show the constraints from iKbl/lKfel and {exl, i-e. 
here and in all the following figures, the region in (p, rj) within the full lines 
is allowed by experiment (see ref. ||14j|). 

Fig. 2b Curves corresponding to constant Q = 0.7,0.85, 1.0, 1.15, 1.28 for 
g = -0.6. 

Fig. 3a Curves corresponding to constant R = 0.025, 0.03, 0.035, 0.04 for 
g = 0.6 and = 0.545. 

Fig. 3b Curves corresponding to constant R = 0.01, 0.03, 0.05, 0.07 for 
g = -0.6 and |#WtP = 0.545. 

Fig. 4a Curves corresponding to constant = 5%, 10%, 15%, 20% for 
g = 0.6, as = 0.2 and K\ = 1.3. 

Fig. 4b Curves corresponding to constant Qcp = 5%, 10%, 15%, 20%, 25%, 
30% for g = -0.6, a, = 0.2 and Ka = 1.3. 

Fig. 5a Curves corresponding to constant N^^ = 2 ■ 10^, 5 ■ 10^, 10^, 5 ■ 10^ 
for g = 0.6, a, = 0.2, Ka = 1.3, A = 0.785 and FpiOj) = -2.6 ■ lO^^^ 

Fig. 5b Curves corresponding to constant N^a = 2 ■ 10^, 5 ■ 10*, 10^, 5 ■ 10^ 
for g = -0.6, a, = 0.2, Ka = 1.3, A = 0.785 and FpiOj) = -2.6 • lO'^. 
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